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Charge-density-wave (CDW) phase transition in a ring-shaped crystals, recently synthesized by 
Tanda et al. [Nature, 417, 397 (2002)], is studied based on a mean-field-approximation of Ginzburg- 
Landau free energy. It is shown that in a ring-shaped crystals CDW undergoes frustration due to the 
curvature (bending) of the ring (geometrical frustration) and, thus, forms a mixed state analogous 
to what a type-II superconductor forms under a magnetic field. We discuss the nature of the phase 
transition in the ring-CDW in relation to recent experiments. 
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I. INTRODUCTION 

Charge-density-wave (CDW) has been attracting in- 
tense interest from solid state physicists, because of its 
unique physical properties as an electronic condensate^. 
Especially analogy to superconductivity has been stimu- 
lating imagination of various people. Frolich paid atten- 
tion to collective transport (sliding) of CDW and inves- 
tigated within it the origin of superconductivity before 
the foundation of BCS theory^. Although it has been 
clarified that the sliding motion of CDW is easily sup- 
pressed by inhomogeneities, realization of the Frolich's 
superconductivity in ultraclean CDW crystals is an in- 
triguing target of experimental trails. 

Similarities and dissimilarities between vortices in su- 
perconductors and dislocations in CDW systems are also 
interesting topics of this field. Analogy between the 
phase slip phenomena via vortex ring in a supercon- 
ducting wire and the dislocation nucleation in a sliding 
CDW has been pointed outi Several people attempted 
to search the counterpart of the mixed state of type-II su- 
perconductors in CDW systems£££ According to Ref. 
la, CDW systems can also form a mixed state analogous to 
that in type-II superconductors when an electric field per- 
pendicular to the chains is applied. This "mixed state" 
of CDW, however, is formed only in the surface region 
of the sample within a depth of Debye screening length, 
which is usually several lattice spacings in most CDW 
materials. Therefore direct experimental observation of 
this state is considered to be very difficult. 

Recently Tanda et al. have succeeded in synthesizing 
a ring crystal of CDW materials 8 , where the conducting 
chains are along the azimuthal direction^S. In a ring 
crystal, CDW system feels frustration due to bending, 
just like a superconductor feels frustration under a mag- 
netic field. Thus we can expect that the CDW forms a 
mixed state as in a type-II superconductor. This fact 
was first pointed out by Nogawa and Nemotoii and the 



dynamics of ring-CDW is investigated in terms of the 
Monte Carlo simulation on a lattice model of CDW. In 
their analysis, ring CDW is mapped to a uniformly frus- 
trated XY model. 



Experimental investigation of the ring-CDW's has 
been carried out by Shimatake et alJ^. They measured 
the single particle relaxation rate in ring- and whisker- 
CDW's using femtosecond spectroscopy and found that 
quasiparticle density in ring CDW's shows a stretched 
exponential decay oc exp{(t/r) - 3 }, whereas it decays ex- 
ponentially in whisker samples. Moreover relaxation time 
in ring CDW's did not show significant divergence at T c , 
which is common in whisker samples. Indication of such 
slow dynamics in ring CDW's is also found in the sim- 
ulations by Nogawa and Nemotoii, who found a power 
law decay of relaxation time in ring crystals. Although 
details of relaxation formula are different, both results 
suggest the existence of slow relaxation processes in ring 
CDW's, which may be due to the dynamics of disloca- 
tions intrinsic to ring CDW's. In this context it is impor- 
tant to study how dislocations intrinsic to ring CDW's is 
distributed in the ordered state. 



In this paper we first derive the Ginzburg-Landau (GL) 
free energy of CDW in a ring geometry. Our results are 
slightly different from uniformly frustrated XY model, or 
superconductors in a magnetic field, especially when the 
radial thickess of the ring is large compared to the radius. 
We then study the mean field behavior of the system 
in order to clarify the transition temperature of CDW 
in a ring geometry. We also comment on the analogy 
to superconductors in a magnetic field, which is useful 
for thin rings (radial thickness -C radius). The relation 
between our results and some experimental results are 
also discussed. 
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II. GINZBURG-LANDAU FREE ENERGY OF A 
RING CDW 

In constructing the free energy describing ring CDW's 
we need detailed knowledge about the crystal structure 
of the ring crystals. We refer to Reffl3l on this point 
and summarize the results. According to this reference, 
crystal defects (edge dislocations) are almost unifromly 
distributed in ring crystals with diameters > 1/xm. This 
means that the chains forming the ring are almost free 
from elastic stress (stretch or contraction) and can be 
described by free CDW chains with different lengths de- 
pending on the radius. There may be some stretch or 
contraction of the order of lattice spacing for each chain 
due to the discrete nature of the lattice, however, this is 
negligible for rings with realistic diameters ~ lpm which 
is much larger than the unit cell of the lattice. 

Then we start from the following GL free energy, 




FIG. 1: Geometry and parameters of our system. The z-axis 
is perpendicular to the sheet. 
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where the first and the second term correspond to intra- 
chain and inter-chain contribution, respectively. The 
number of the chains is denoted by N. Summation over k 
corresponds to the stacking of the layer along the axis of 
the ring, with N z being the number of the layers. Since 
we assume that the system is uniform along the axis, k- 
summation is suppressed hereafter. Here is the free 
energy of the J-th chain given by 
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where Rj is the radius of the j-th chain given by Rj = 
R + d (i — "2") (R ^ s the radius at the center of the width 
and d is the inter-chain spacing), £o is the coherence 
length of CDW at zero temperature T = 0, and t is the 
reduced temperature T/T c with T c being the transition 
temperature in the bulk. ipj( x ) is the order parameter 
on the j-th chain. 

The inter-chain coupling term is constructed in the 
following way: We write the order parameter as ?pj = 
|^,.| e »0j(z) _ Then the local charge density modula- 
tion due to CDW can be written as 6p [ ° cal (x) oc 
\ipj{x)\cos{2kpx + (f>j(x)}, and the inter-chain coupling 
is given in the following form, 
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where Jo is a constant. CDW usually prefers 7r-phase 
shift between neighboring chains due to, for example, 
Coulomb repulsioni. We have beforehand taken it into 



account in the definition of the order parameter, and the 
minus sign of Eq. ([3]) is due to this. Noting the equality, 

cos{2kpRj9 + 4>j(x)} cos{2fcpi? J+ i# + <j)j +1 (x)} 

= I cos{2k F (Rj + R j+l )6 + 4j(x) + 4> 3+ i{x)} 



+ i cos{2k F d9 + 4>j{x) - (f> j+ i(x)}, 



(4) 



where d = Rj+i — Rj, and neglecting the first term of 
the right hand side as the negligible fast oscillation, we 
obtain, 
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In the last line, we added some terms proportional to 
\ipj\ 2 , in order to subtract the effect of the renormaliza- 
tion of T c . 

Next we gauge transform the order parameter as 
ijjj{x) — > ipj(x)e z2kF J and obtain the resulting free en- 
ergy as, 
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Since we are interested in the behavior of CDW near 
T c , we have disregarded the effects of electrostatic screen- 
ing. In order words, the scalar potential does not appear 
in our equations. This is an over simplification when we 
go to lower temperatures, where CDW order is well de- 
veloped and charge modulation is significant. This point 
will be discussed further in the end of this paper. 

The boundary condition for ipj (x) in the ring-geometry 
is given by 



smaller lengths. This can make the thermodynamic prop- 
erties of ring CDW's different from those in superconduc- 
tors in a magnetic field especially when the ring is thick 
W ~ R. This point should be carefully studied in further 
analysis. 



III. CASE OF A THIN RING W < R: A 
SUPERCONDUCTING ANALOGY 



d x ipj{0) = d x ipj{Rj). 



(8) 
(9) 



These conditions impose 4>j{Rj) — <fij(0) = —^nkpRj + 
2ttN, where N is an integer. This gives rise to the phe- 
nomena similar to the Little-Parks effect in supercon- 
ducting rings. For systems much larger than the Fermi 
wave length, Rj ^ l/kp, this effect is very weak and 
negligible. This applies for most ring CDW systems syn- 
thesized until now. 

Next we consider a change of coordinate (from chain 
to chain) x — > x', so that Rj becomes R for all the chains 
(Namely, x/Rj = x'/R). Here defining x{Rj) = — 
R)k,F / R, we obtain 
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Although thick rings (W ~ R) are a little different 
from superconductors in a magnetic field, as we have seen 
above, the analogy to superconductors is still useful for 
thin rings (W <C R). For thin rings, the free energy is 
simplified to 
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where 7 is defined by Jo/4 = l 2 Q/d 2 . This is actually 
the continuum version of the lattice model previously in- 
troduced by Nogawa and Nemoto^. 

We further change the coordinates as y' = y/j so that 
the system is isotropic in both x- and y-direction and 
obtain 



The inter-chain term can be rewritten as 
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Now we apply continuum approximation to the radial 
direction (y-direction). Putting W — Nd, we obtain 
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It is worthwhile to point out that the free energy F = 
pc _|_ pmt jjg^g j-affogj- complicated y dependences due 
to the bending. It is not simply the type of the frus- 
trated XY model, or a superconductor in a magnetic 
field. Especially, the factor (y/R) appearing before x'- 
integral comes from the fact that inner chains have 
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with 



e x being the unit vector in x-direction. 

At this stage, we see a nice analogy to the free energy 
of a superconductor in a magnetic field. The magnetic 
field H corresponds to the curvature of the ring 1 /R as 
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From this analogy we can deduce the transition temper- 
ature of CDW in a ring crystal. We note that a ring 
CDW corresponds to a type-II superconductor, not type- 
I, since there is no mechanism to screen the "curvature" 
(or bending) of the ring in the CDW system. In other 
words, there is no Meissner phase. 

Since the transition temperature of a type-II superd- 
conductor under a magnetic field H is given by 



T c2 (H) = T c (0) 1 
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where H C 2(0) = (j>o / '(27t£q) is the upper critical field at 
T = (Note that this formula is valid if H < H c2 (0). In 
that sense, H C 2(Q) is the critical field at T = obtained 
from the linear extrapolation from T ~ r c .), the CDW 
transition temperature of a ring crystal of radius R is 
given by 



(18) 



where i? C 2 = 2&f7£o and T c (oo) is the transition tem- 
perature of the whisker sample. 

Here we apply our results to actual NbSe3 ring 
crystal. We put 2kp — 4.3nm _1 , £o — 6.1nm and 
7 ~ 1/10 ~ 1/100 for NbSe 3 , and then we obtain 
R C 2 ~ 160 x 7nm = 1.6 ~ 16nm. 



IV. NUMERICAL STUDY OF LINEARIZED 
GINZBURG-LANDAU EQUATION 

We study the linearized GL equation for the cases of 
arbitrary R and W and determine the CDW transition 
temperature. We start by disregarding the quartic term 
of the order parameter in Eqs. (fT0|) . By introducing the 
Fourier transformation, 
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We put ^-(n) = y/Rj/Ripjfo), and the above equation 
becomes 
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The integer n indexes the mode of the order parameter. 
All the eigenvalues of Mj&(n) are positive and we de- 
note the smallest one by A m i n (n). The mode indexed by 
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text) 
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0.5aRg (see 



n becomes unstable when A m i n («) + t — 1 < is satis- 
fied. Below this temperature ^j(n) (or ipj(n)) tends to 
have a nonzero expectation value and the nonlinear term 
of GL free energy becomes relevant. The mode which 
becomes unstable at the highest temperature dominates 
the real CDW transition. Therefore, by diagonalizing 
Mjk(n) and comparing the eigenvalues of all possible n's, 
we can estimate the CDW transition temperature. 

In Figs. [2] and [3j we have shown the log-log (natural 
logarithm) plot of the shift of the transition temperatures 
of several systems as a function of the radius of outermost 
chain R OM t = R + W/2, (AT C is the decrease of the 
transition temperature from the bulk T c .) Fig. [3 is for 
thick samples where W ~ 2R (W = 2(R - d) to be 
more precise). In this case, the ring looks almost like a 
disk. Fig. [3] is for thin samples where we put W = R/5. 
(a) and (b) in each figure are data for different coherence 
length £o- (a) is for Co = 5d and (b) is for £ = 20c?. From 
these figures, we can see that amount of AT C is inversely 
proportional to R ou t for samples with larger diameters. 

We write T c (R out ) — T c (l — R c /R out ) and estimate R c 
from the numerical data. The results are shown in the 
table [H R c is obtained from the fitting in the figures. If 
the fitting line is y — ax + where a ~ —1, we obtain 
R c = e 13 . It can be seen that R c is not close to R C 2 which 
is estimated from Eq. (fT5)l . We found that R c is rather 
close to R C 3 which corresponds to the third critical field 
H C 3 of type-II superconductors. According to Ref. 0, 
R c2 and R c3 are related by R^ = 0.59R% H c3 of a 
type -II superconductor gives the critical field of surface 
nucleation of superconducting order in a magnetic field. 
Therefore R c in a ring CDW may also correspond to the 
2 same phenomenon. This is confirmed by examining the 
order parameter which develops at the transition tem- 
perature. In Fig. 21 we have shown the order parameter 
at the critical temperature in case of Fig. [2] (b). Even 
though the radius is changed, CDW order starts devel- 
oping always from the outer edge of the sample and it 
tells us that the CDW transition is of surface nucleation 
type. 

From the above argument we obtain a rough picture 
of CDW transition in the ring crystals: At T — T c3 = 
T c (l — R c3 /R), CDW order starts nucleation at the outer 
surface of the sample and, at T = T c2 = T c (l — R c2 /R), 
CDW order prevails in the entire sample. 
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FIG. 2: Log-log plot of the outermost radius R ou t vs. the 
shift of the transition temperature in thick (W ~ 27?) rings 
(or cylinders), (a) £o = 5d, (b) £o = 20d. Lines are the fitting 
line to the linear part of the data. Rout is in the unit of d. 
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FIG. 3: Log-log plot of the outermost radius Rout vs. the 
shift of the transition temperature in thin (W = R/5) rings, 
(a) £o = 5d, (b) £o = 20d. Lines are the fitting line to the 
linear part of the data. R ou t is in the unit of d. 



V. DISCUSSION 

Let us compare our theory with the experimental re- 
sults obtained by Tsuneta et alJ^ 

Tsuneta et al. found that the transition tempera- 
ture estimated from the anomalies in resistivity does not 
change significantly between whisker and ring crystals, 
although significant decrease is observed in "figure-eight" 
(twisted ring) cryslats. As for the ring crystals, their 
data are in good agreement with our result, since the ra- 
diuses of the ring crystals used in the experiment may 
be mostly above 10/im and then, according to our the- 
ory, the suppression of transition temperature is about 
T c R c2 /R - 0.001 x T c - 0.1K. This is smaller than the 
sample-to-sample fluctuations. In contrast to this, in 
figure-eight crystals, the transition temperature is sup- 
pressed about IK. This significant suppression compared 
to rings is worth noting. Although detailed effects of 
the twisting is not clear at this stage, we consider that 
such a large suppression of T c due to geometrical frustra- 
tion is unexpected from the naive intuition based on the 
present analysis. The explanation of this phenomenon 
needs further studies, both from theoretical and experi- 
mental sides. 



Here we comment on the lower temperature side. The 
free energy, Eqs. (|10p and pip , is valid only in very vicin- 
ity of T c , since the effect of electrostatic screening (Debye 
screening) is neglected. The dislocations in CDW give 
rise to the polarization of charge arising from d x (f>j(x), 
and thus undergoes the electrostatic screening due to the 
quasiparticles and the condensate (phasons). Since the 
electrostatic screening is strong in most CDW materials, 
this neglect is serious. A free energy which includes the 
screening effect is given, for example, in Ref. 0, accord- 
ing to which, the interaction between dislocations is not 
the type of anisotropic GL theory but shows complicated 
dependence on direction and scale. In order to clarify the 
CDW order in ring crystals at much lower temperature 
than T c , a study based on such accurate free energy is 
indispensable, which may need further numerical studies. 
We leave this for future study. 

It should also be noted that in realistic CDW materials 
the effects of thermal fluctuation are essential. As is well 
known, the real transition temperature of CDW is much 
lower than the mean-field value which is estimated from 
the CDW energy gap. This is considered to be due to the 
renormalization of T c by the thermal fluctuations. In or- 
der to obtain more reliable results, we need to go beyond 
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FIG. 4: The amplitude of the order parameter at the CDW 
transition temperature for samples with several different ra- 
diuses. Each peak corresponds to solutions for different ra- 
dius. The peak is always located at the outer edge of the 
sample. The inset is the magnification of one of the data. 
Bold line shows the outer edge. 



the mean field approximation, For example, numerical 
simulations as performed by Nogawa and Nemoto, given 
in Ref. [HI, are hopeful. Further studies in such direction 
are desired. 



VI. SUMMARY 

We have derived the free energy which describes the 
CDW order in ring crystals and, based on it, calcu- 
lated the mean-field transition temperature. We have 
discussed several physical properties of ring CDW based 
on the obtained free energy. 
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